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Scaling laws that describe the structure of river networks are shown to follow from three simple 
assumptions. These assumptions are: (1) river networks are structurally self-similar, (2) single chan- 
nels are self-affine, and (3) overland flow into channels occurs over a characteristic distance (drainage 
density is uniform). We obtain a complete set of scaling relations connecting the exponents of these 
scaling laws and find that only two of these exponents are independent. We further demonstrate 
that the two predominant descriptions of network structure (Tokunaga's law and Horton's laws) are 
equivalent in the case of landscapes with uniform drainage density. The results are tested with data 
from both real landscapes and a special class of random networks. 

PACS numbers: 92.40.Fb, 92.40.Gc, 68.70.+W, 64.60. Ht 



I. INTRODUCTION 

If it is true that scaling laws abound in nature jl), 
then river networks stand as a superb epitome of this 
phenomenon. For over half a century, researchers have 
uncovered numerous power laws and scaling behaviors 
in the mathematical description of river networks 
% [§ |, % 0' I' D- Tnese scaling laws, which are 
usually parameterized by exponents or ratios of funda- 
mental quantities, have been used to validate scores of 
numerical and theoretical models of landscape evolu- 
tion 0, 0, g [D§ [fjj @ H [f| g, j§ H 
and have even been invoked as evidence of self-organized 
criticality [^3], However, despite this widespread us- 
age, there is as yet no fundamental understanding of the 
origin of scaling laws in river networks. 

It is the principal aim of this paper to bring together 
a large family of these scaling laws within a simple, logi- 
cal framework. In particular, we demonstrate that from 
a base of three assumptions regarding network geome- 
try, all scaling laws involving planform quantities may 
be obtained. The worth of these consequent scaling laws 
is then seen to rest squarely upon the shoulders of the 
structural assumptions themselves. We also simplify the 
relations between the derived laws, demonstrating that 
only two scaling exponents are independent. 

The paper is composed in the following manner. We 
first present preliminary definitions of network quantities 
and a list of empirically observed scaling laws. Our as- 
sumptions will next be fully stated along with evidence 
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for their validity. Several sections will then detail the 
derivations of the various scaling laws, being a combina- 
tion of both new insights of our own as well as previous 
results. Progressing in a systematic way from our as- 
sumptions, we will also be required to amend several in- 
consistencies persistent in other analyses. The theory will 
be tested with comparisons to data taken from real land- 
scapes and Scheidegger's random network model pa, p6[. 



II. THE ORDERING OF STREAMS 

A basic tool used in the analysis of river networks 
is the device of stream ordering. A stream ordering is 
any scheme that attaches levels of significance to streams 
throughout a basin. Most orderings identify the small- 
est tributaries as lowest order streams and the main or 
'trunk' stream as being of highest order with the inter- 
mediary 'stream segments' spanning this range in some 
systematic fashion. Stream orderings allow for logical 
comparisons between different parts of a network and 
provide a basic language for the description of network 
structure. 

Here, we build our theory using the most common or- 
dering scheme, one that was first introduced by Horton in 
his seminal work on erosion j2|. Strahler later improved 
this method p7| and the resulting technique is commonly 
referred to as Horton-Strahler stream ordering p~3| . The 
most natural description of this stream ordering, due to 
Melton p8| , is based on an iterative pruning of a tree rep- 
resenting a network as shown in Figure ^. All source (or 
external) streams are pared away from the tree, these be- 
ing defined as the network's first order 'stream segments'. 
A new tree is thus created along with a new collection 
of source streams and these are precisely the second or- 
der stream segments of the original network. The prun- 
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ing and order identification continues in like fashion until 
only the trunk stream segment of the river network is left. 
The overall order of the basin itself is identified with the 
highest stream order present. 

The usual and equivalent description details how 
stream orders change at junctions J23|. When a stream 
segment of order lo\ merges with a stream segment of or- 
der 0J2 , the outgoing stream will have an order of u given 
by 

t l ) = max(ui,w 2 )+i 1 , lJ2 (1) 

where 5 is the Kronecker delta. In other words, stream 
order only increases when two stream segments of the 
same order come together and, otherwise, the highest 
order is maintained by the outflowing stream. 



III. PLANFORM NETWORK QUANTITIES 
AND SCALING LAWS 

The results of this paper pertain to networks as viewed 
in planform. As such, any effects involving relief, the 
vertical dimension, are ignored. Nevertheless, we show 
that a coherent theory of planform quantities may still 
be obtained. This section defines the relevant quantities 
and their various permutations along with scaling laws 
observed to hold between them. The descriptions of these 
laws will be short and more detail will be provided in later 
sections. 

The two essential features in river networks are basins 
and the streams that drain them. The two basic plan- 
form quantities associated with these are drainage area 
and stream length. An understanding of the distribu- 
tion of these quantities is of fundamental importance in 
geomorphology. Drainage area, for example, serves as a 
measure of average discharge of a basin while its relation- 
ship with the length of the main stream gives a sense of 
how basins are shaped. 



A. General network quantities 

Figure || shows a typical drainage basin. The basin 
features are a, the area, Z, the length of the main stream, 
and L\\ and Lj_, the overall dimensions. The main (or 
trunk) stream is the dominant stream of the network — 
it is traced out by moving all the way upstream from 
the outlet to the start of a source stream by choosing 
at each junction (or fork) the incoming stream with the 
largest drainage area. This is not to be confused with 
stream segment length which only makes sense in the 
context of stream ordering. We will usually write L for 
L\\. Note that any point on a network has its own basin 
and associated main stream. The sub-basin in Figure ^ 
illustrates this and has its own primed versions of a, I, 
Lu and Lj_- The scaling laws usually involve compar- 
isons between basins of varying size. These basins must 




(b) (c) 

FIG. 1: Horton-Strahler stream ordering, (a) shows the ba- 
sic network, (b) is created by removing all source streams 
from the network in (a) , these same streams being denoted as 
first order 'stream segments'. The new source streams in the 
pruned network of (b) are labelled as second order stream seg- 
ments and are themselves removed to give (c), a third order 
stream segment. 

be from the same landscape and may or may not be con- 
tained within each other. 

Several scaling laws connect these quantities. One of 
the most well known is Hack's law ||. Hack's law states 
that I scales with a as 

I ~ a h (2) 

where h is often referred to as Hack's exponent. The 
important feature of Hack's law is that h ^ 1/2. In 
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particular, it has been observed that for a reasonable 
span of basin sizes that 0.57 < h < 0.60 [§, |, || |o). 
The actual range of this scaling is an unresolved issue 
with some studies demonstrating that very large basins 
exhibit the more expected scaling of h = 1/2 |n], ^2|, ^3) . 
We simply show later that while the assumptions of this 
paper hold so too does Hack's law. 

Further comparisons of drainage basins of different 
sizes yield scaling in terms of L(= Ln), the overall basin 
length. Area, main stream length, and basin width are 
all observed to scale with L @ § § §, ||, 

a~L D , l~L d , L± ~ L H . (3) 

Turning our attention to the entire landscape, it is also 
observed that histograms of stream lengths and basin 
areas reveal power law distributions (9[ ^3) : 

P(a) ~ a- T and P(l) ~ r 7 . (4) 

There are any number of other definable quantities and 
we will limit ourselves to a few that are closely related 
to each other. We write A for the average distance from 
a point on the network to the outlet of a basin (along 
streams) and A for the unnormalized total of these dis- 
tances. A minor variation of these are A and A, where 
only distances from junctions in the network to the outlet 
are included in the averages. 

The scaling law involving these particular quantities is 
Langbein's law H which states that 

A ~ a' 3 . (5) 
Similarly, we have A - LP, A - aP and A ~ L^, [§. 




FIG. 2: A planform view of an example basin. The main 
defining parameters of a basin are a, the drainage area, I, 
the length of the main stream, and L« and L±, the overall 
Euclidean dimensions. The sub-basin with primed quantities 
demonstrates that a basin exists at every point in a network. 



where Rt is a fixed constant for a given network. Thus, 
all of Tokunaga's ratios may be specified by two funda- 
mental parameters T\ and Rt'- 



B. Network quantities associated with stream 
ordering 



T^T^Rt)" 



(7) 



With the introduction of stream ordering, a whole new 
collection of network quantities appear. Here, we present 
the most important ones and discuss them in the context 
of what we identify as the principal structural laws of 
river networks: Tokunaga's law and Horton's laws. 



1. Tokunaga's law 

Tokunaga's law concerns the set of ratios, {T^/}, 
first introduced by Tokunaga |||, [56|, [57| |3^, p9[ . 
These 'Tokunaga ratios' represent the average number of 
streams of order to' flowing into a stream of order to as side 
tributaries. In the case of what we will call a 'structurally 
self-similar network', we have that T w%w i = T^^^i = T v 
where v = lo — u>' since quantities involving comparisons 
between features at different scales should only depend 
on the relative separation of those scales. These T v , in 
turn, are observed to be dependent such that pjl , 



T v +i/T u = Rt 



(6) 



We refer to this last identity as Tokunaga's law. 

The network parameter T\ is the average number of 
major side tributaries per stream segment. So for a col- 
lection of stream segments of order w, there will be on 
average T\ side tributaries of order u> — 1 for each stream 
segment. The second network parameter Rt describes 
how numbers of side tributaries of successively lower or- 
ders increase, again, on average. As an example, consider 
that the network in Figure l] is part of a much larger net- 
work for which T± = 2 and Rt — 4. Figure [l] (b) shows 
that the third order stream segment has two major side 
tributaries of second order which fits exactly with T± = 2 
(Note that the two second order stream segments that 
come together to create the third order stream segment 
are not side tributaries). Figure [l] (a) further shows nine 
first order tributaries, slightly above the average eight 
suggested by T-2 = T\R T = 8. Finally, again referring to 
Figure |l] (a), there are 9/4 = 2.25 first order tributaries 
for each second order stream segment, not far from the 
expected number Ti = 2. 
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2. Morton's laws 

Horton introduced several important measurements for 
networks in conjunction with his stream ordering ||] . The 
first is the bifurcation ratio, R n . This is the ratio of the 
number n u of streams of order uj to the number 
of streams of order uj + 1 and is, moreover, observed to 
be independent of ui over a large range. There is next 
the stream length ratio, Re = where £ u is the 

average length of stream segments of order u>. These 
lengths only exist within the context of stream ordering. 
In contrast to these are the main stream le ngths, which 
we have denoted by I and described in section III A . Main 
stream lengths are defined regardless of stream ordering 
and, as such, are a more natural quantity. Note that 
stream ordering gives rise to a discrete set of basins, one 
for each junction in the network. We therefore also have 
a set of basin areas and main stream lengths defined at 
each junction. Taking averages over basins of the same 
order we have a w and l u to add to the previously defined 
l u and n u . 

The connection between the two measures of stream 
length is an important, if simple, exercise fllpf . Assuming 
£uj+i = Riluj holds for all uj, one has 



Law: 



Name or description: 



T v = Ty(Rt) Tokunaga's law 

I ~ L d self-affinity of single channels 

n^+i/n^ = R n Horton's law of stream numbers 

= Re Horton's law of stream segment lengths 

lu+i/lu ~ Re Horton's law of main stream lengths 

a u +i/fl u = R a Horton's law of stream areas 



I - 


a h 


Hack's law 


a ~ 


L D 


scaling of basin areas 




^L H 


scaling of basin widths 


P(a) 


~ a~ T 


probability of basin areas 


P(l) 


~ r 1 


probability of stream lengths 


A r- 




Langbein's law 


A - 




variation of Langbein's law 


A - 


-a* 


as above 


A - 




as above 



TABLE I: A general list of scaling laws for river networks. All 
laws and quantities are defined in section HI. The principal 
finding of this paper is that these scaling laws follow from the 
first two relations, Tokunaga's law (structural self-similarity) 
and the self-affinity of single channels, and t he ass umption of 
uniform drainage density (defined in section IV C ) . 



2=1 



[Re 



1 



Re -l 



(8) 



C. Scheidegger's random networks 



where h = l\ has been used. Since typically Rt > 2 |glj, 
— > Re rapidly. For uj = 4 and Re — 2, the error 
is only three per cent. On the other hand, starting with 
the assumption that main stream lengths satisfy Horton's 
law of stream lengths for all uj implies that the same is 
true for stream segments. 

Thus, for most calculations, Horton's law of stream 
lengths may involve either stream segments or main 
streams and, for convenience, we will assume that the 
law is fully satisfied by the former. Furthermore, this 
small calculation suggests that studies involving only 
third- or fourth-order networks cannot be presumed to 
have reached asymptotic regimes of scaling laws. We will 
return to this point throughout the paper. 

Schumm p2| is attributed with the concrete introduc- 
tion of a third and final law that was also suggested by 
Horton. This last ratio is for drainage areas and st ates 
that R a — au+i/a^- We will later show in section VII 
that our assumptions lead to the result that R a = R n . At 
this stage, however, we write Horton's laws as the three 
statements 



To end this introductory section, we detail some of the 
features of the random network model of Scheidegger |2f| 
Although originally defined without reference to a 
real surface, Scheidegger networks may be obtained from 
a completely uncorrelated landscape as follows. Assign 
a random height between and 1 at every point on a 
triangular lattice and then tilt the lattice so that no local 
minima (lakes) remain. Scheidegger networks are then 
traced out by following paths of steepest descent. 

Surprisingly, these networks still exhibit all of the scal- 
ing laws observed in real networks. It thus provides an 
important point in 'network space' and accordingly, also 
provides an elementary test for any theory of scaling laws. 
Exact analytical results for various exponents are known 
due to the work of Takayasu et al. on the aggregation of 
particles with injection Q |||, ||, |46[ Q . While 
there are no analytic results for the Tokunaga ratio T\ or 
the Horton ratios R n and Re, our own simulations show 
that these stream order laws are strictly obeyed. Table || 
lists the relevant exponents and their values for the Schei- 
degger model along with those found in real networks. 



— R n , — =~— — Re j snd _^ — R a . 



(9) 



IV. ASSUMPTIONS 



A summary of all of the scaling laws presented in this 
section is provided in Table Q. Empirically observed val- 
ues for the relevant exponents and ratios are presented 
in Table ||. 



We start from three basic assumptions about the struc- 
ture of river networks: structural self-similarity, self- 
affinity of individual streams and uniformity of drainage 
density. We define these assumptions and their relevant 
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Quantity: Scheidegger: Real networks: 



Rn 
Ra 

Rt 

d 
D 
h 

T 

7 
¥> 
H 

P 
<P 



5.20 ± .05 
5.20 ± .05 
3.00 ± .05 
1.30 ±.05 
1 

3/2 
2/3 
4/3 
3/2 

1 
1/2 
5/3 

1 



3.0-5.0 4£ 
3.0-6.0 \~ 
1.5-3.0 1 49 
1.0-1.5 p 
1.1 ±0.01 
1.8 ±0.1 
0.57-0.60 
1.43 ±0.02 
1.8 ±0.1 
1.05 ±0.01 
0.75-0.80 
1.56 | 
1.05 ±0.01 



TABLE II: Ratios and scaling exponents for Scheidegger's 
random network model and real networks. For Scheidegger's 
model, exact values are known due to the work of Takayasu et 

ai. @ 0, m n || n and approximate results are taken 

from our own simulations. For real networks, the references 
given are generally the most recent and further appropriate 
references may be found within them and also in section M . 



parameters and then discuss their mutual consistency. 
We end with a discussion of the correspondence between 
the laws of Tokunaga and Horton. It should be stressed 
that while we make a case for each assumption there is 
also considerable proof to ponder in the pudding that 
these ingredients create. 



A. Structural self-similarity 

Our first assumption is that networks are structurally 
self-similar. It has been observed that river networks 
exhibit self-similarity over a large range of scales ||, 
|23| . Naturally, the physical range of this self-similarity 
is restricted to lie between two scales. The large scale 
cutoff is the overall size of the landscape and the small 
scale cutoff is of the order of the characteristic separation 
of channels |5(| . 

In order to quantify this phenomenon, we look to laws 
of network structure such as Tokunaga's law and Hor- 
ton's laws of stream number and length. We demon- 
strate in the following section that these descriptions are 
mutually consistent within the context of our third as- 
sumption, uniformity of drainage density. Thus, we may 
assume a network where both Tokunaga's and Horton's 
laws hold. For convenience, we write these laws as if they 
hold for all orders down to the first order. Any actual 
deviations from these laws for low orders will not affect 
the results since we are interested in how laws behave for 
increasing stream order. 



B. Self-affinity of individual streams 

Our second assumption is that individual streams are 
self-affine curves possessing a dimension d > 1, as intro- 
duced in equation (0) . Empirical support for this premise 
is to be found in m, || , [| ^3|, [|l[] . In reality, this is at 
best a weak fractality with measurements generally find- 
ing d to be around 1.1 ||. We assume d to be constant 
throughout a given network, true for each stream inde- 
pendent of order. 

In general, it is most reasonable to consider this in the 
sense of a growing fractal: stream length I will grow like 
L d where L is the overall length of a box containing a 
portion of a stream. So, rather than examine one fixed 
section of a stream, we take larger and larger pieces of it. 
Moreover, this is the most reasonable method for actually 
measuring d for a real network. 



C. Uniform drainage density 

Our third and final assumption is that drainage den- 
sity is uniform throughout a network. For a given basin, 
the drainage density, p, is a measure of the average 
area drained per unit length of stream by overland flow 
(i.e., excluding contributions from tributary streams). Its 
usual form is that given by Horton B : 



P : 



(10) 



where, for a given basin, ^ t represents the summed to- 
tal length of all stream segments of all orders and a is 
the drainage area. More generally, one can in the same 
way measure a local drainage density for any connected 
sections of a network within a landscape. Such sections 
should cover a region at least l\ in diameter, the typical 
length of a first order stream. Drainage density being 
uniform means that the variation of this local drainage 
density is negligible. There is good support in the liter- 
ature for the uniformity of drainage density in real net- 
works H [|| H| i§) while there are some sugges- 
tions that it may vary slightly with order |3^] . 

Uniform drainage density may also be interpreted as 
the observation that the average distance between chan- 
nels is roughly constant throughout a landscape p3[ , 
an estimate of this distance being simply 1/ p. This is due 
to the fact that there is a finite limit to the channelization 
of a landscape determined by a combination of soil prop- 
erties, climate and so on. Implicit in this assumption is 
that the channel network has reached its maximum ex- 
tension into a landscape JE^, [57}. Indeed, In the bold 
words of Glock J57|, we are considering river networks 
at the "time of completed territorial conquest." Further- 
more, Shreve Q notes that drainage density would be 
uniform in a "mature topography developed in a homo- 
geneous environment." 

Importantly, our third assumption connects the plan- 
form description to the surface within which the network 
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lies. Computationally, the uniformity of drainage density 
allows for the use of the length of a stream as a proxy for 
drainage area |Efj|. Further, the average distance between 
streams being roughly constant implies that, on average, 
tributaries are spaced evenly along a stream. 



V. TOKUNAGA'S LAW AND HORTON'S LAWS 
ARE EQUIVALENT 

This section demonstrates an equivalence between 
Tokunaga's law and Horton's two laws of stream num- 
ber and stream length in the case of a landscape with 
uniform drainage density. 



A. From Tokunaga's law to Horton's laws 

Tokunaga has shown that Horton's law for stream 
numbers follows from Tokunaga's law (given in equa- 
tion (Q)) |}6], [38|. This follows from the observation that 
n u , the number of streams of order w, in a basin of order 

may be expressed as 



of lengths of stream between tributaries is then simply 
one more in number. 

Using Tokunaga's law (equation (Q)) we find that 



L+i/L = r t (i + o(r t )-") 



(14) 



obtaining Horton's stream length ratio with the simple 
identification: 



Ri = Rt 



(15) 



and we will use Re in place of Rt throughout the rest of 
the paper. As already noted we will see that R a = R n 
for landscapes where drainage density is uniform. This 
redundancy means that there are only two independent 
Horton ratios, Ri and R n , which sits well with the two 
independent quantities required for Tokunaga's law, T\ 
and Rt- Presupposing this result, we can invert equa- 
tions ( |l2] ) and (|l5|) to obtain Tokunaga's parameters from 
the two independent Horton ratios: 



Rt 
Ti 



Re 

Rn 



Re ~ 2 + 2Re/R ri 



(16) 
(17) 



2n 



(11) 



B. From Horton's laws to Tokunaga's law 



The 2n w +i accounts for the fact that each order lo + 1 
stream is initiated by the confluence of two streams of 
order lo. Presuming Tokunaga's law, a simple analysis of 
equation ( pd| ) shows that in the limit of large f2, the ratio 
n u /n u+ i does indeed approach a constant. This leads to 
an expression for the Horton ratio R n in terms of the 
two Tokunaga parameters T\ and Rt (first obtained by 
Tokunaga in |3q|): 



2R n 



(2 + R T + T 1 )+[(2 + R T + T 1 ) 2 -8R T ] 1/2 ■ (12) 



Tokunaga's work has been recently generalized by 
Peckham who deduces links to the other Horton ratios 
Re and R a [p8| . In contrast to the purely algebraic cal- 
culation of R n , these results require the step of equat- 
ing topological properties to metric basin quantities. In 
determining Re, Peckham uses the number of side trib- 
utaries to a stream as an estimate of stream segment 
length. This is based on the assumption t hat t ributaries 
are evenly spaced. As discussed in section [V C , this even 



spacing of tributaries follows for networks with uniform 
drainage density. Therefore, we may write, after Peck- 
ham, that 



L oc 1 + 



E 



T 



(13) 



where the dimension of length absent on the right-hand 
side is carried by an appropriate constant of proportion- 
ality. This sum is simply the total number of tributaries 
that, on average, enter a stream of order lo. The number 



We now provide an heuristic argument to show that 
Tokunaga's law in the form of equation (0) follows from 
Horton's laws of stream number and length and uniform 
drainage density. Note that even though we have shown 
in equations ([l2|), (|l5|), and ([l7]) that the parameters of 
Tokunaga's law and those of Horton's laws may be ob- 
tained from each other, it is not a priori clear that this 
result would be true. Indeed, Tokunaga's law contains 
more direct information about network structure than 
Horton's laws and it is the additional constraint of uni- 
form drainage density that provides the key. 

Consider a stream of order to along with its side trib- 
utaries of order lo' = 1 through co' = lo — 1, the numbers 
of which are given by the usual T v where v = lo — lo' 
(see Figure ||). Since the presumed adherence to Hor- 
ton's laws implies that a network is self-similar we need 
only consider the form of the T v and not the more general 
T u ' jU . Now, again since networks are self-similar, a typi- 
cal stream of order lo + 1 can be obtained by scaling up 
the picture of this order lo stream. As per Horton's law 
of stream lengths, this is done by increasing the length 
of each stream by a factor of Re (Figure || (a) becomes 
Figure 1(b)). 

However, since order lo' streams become lo' + 1 streams 
in this rescaling, the picture in Figure |^ (b) is missing 
first order streams. Also, the average distance between 
tributaries has grown by a factor of Re- Therefore, to 
retain the same drainage density, an extra (Re — 1) first 
order streams must be added for each link (one more than 
the number of tributaries) along this new order lo + 1 
stream (Figure (c)). Since the number of first order 
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(a) 



(b) 



(c) 




FIG. 3: An example rescaling of a basin to demonstrate 
how Tokunaga's law follows from Horton's laws and uniform 
drainage density. In the first step from (a) to (b), the streams 
of the small network are rescaled in length by a factor of 
Re. The second step from (b) to (c) demonstrates that for 
drainage density to remain constant and uniform, a sufficient 
number of first order tributaries must be added. 



streams is now given by T w+ i we have 

T u+1 = {R e - 1) r£ T„ + lj. (18) 

It may be simply checked that this equation is satisfied, 
for large u>, by Tokunaga ratios given by equation (0). 
Thus, Horton's laws of stream number and stream length 
and the uniform drainage density are seen to imply Toku- 
naga's law. 

In general, Horton's ratios rather than the parameters 
of Tokunaga's law will be the most useful parameters 
in what follows. In particular, we will see that the two 
independent quantities R n and Re will be needed only 
in the form lni?„/lni?f. All other exponents will be 
expressible as algebraic combinations of In R n / In Re and 
d, the fractal dimension of an individual stream. 

Furthermore, example (or modal) values for the pa- 
rameters of Horton and Tokunaga are ||(]. fll]] 

Ti = 1, R T = Re = 2, and R n = 4. (19) 

The parameters have been chosen so as to satisfy the 
inversion relations of equation (|l^) . As shown in Table 0, 
real networks provide some variation around these modal 
values. These will be used as rough checks of accuracy 
throughout the rest of the paper. 



VI. HACK'S LAW 

One of the most intriguing scalings found in river net- 
works is Hack's law Q which relates main stream length 
to basin area as I ~ a h . This equation has been empir- 
ically shown to hold true for a large range of drainage 
basin sizes on many field sites p3| . The salient feature 
is that for smaller basins |3(|, h is typically found to be 
in the range (0.56,0.60), whereas 0.5 would be expected 
from simple dimensional analysis [p3| . 

It should be emphasized that Hack's law is only true on 
average as are, for that matter, Tokunaga's law and Hor- 
ton's laws. An extension of Hack's law to a more natural 
statistical description of the connection between stream 
lengths and drainage areas was suggested by Maritan ct 
al. p] with some further developments to be found in [p8[. 



A. Horton's other law of stream numbers 

In order to obtain Hack's law, we will use the unifor- 
mity of drainage density to estimate the area of an order 
n basin by calculating the total length of streams within 
the same basin. So we simply need the typical length 
and number of each stream order present. Taking the 
length of a source stream, ti, to be the finest resolution 
of the network and the basic unit of length, the length 
of a stream segment of order ui is = {ReY~ x l\. How- 
ever, in finding the frequency of such streams we find 
that some care must be taken for the following reasons. 

Horton's law of stream numbers is potentially mislead- 
ing in that it suggests, at first glance, that within a basin 
of order u> there should be one stream of order u), R n 
streams of order u> — 1, R\ streams of order uj — 2 and 
so on. Indeed, many calculations involving Horton's laws 
use this assumption |?], |2^, |5(| [5!| . 

But Horton's R n actually provides the ratio of the 
number of streams of consecutive orders as totalled for 
a whole basin. To illustrate this fact, consider streams 
of order tu and uj + 1 within a basin of order Q 3> u). 
As Tokunaga's law makes clear, streams of order uj are 
not all found within sub-basins of order u> + 1. Indeed, 
a certain number of order uj streams will be tributaries 
to streams of order greater than uj + 1 (see the example 
network of Figure [j] (a)). Tokunaga's law shows that we 
should in fact expect T\ + 2 rather than than R n streams 
of order uj entering into a stream of order u + 1. For 
the typical values T\ = 1 and R n = 4 in (|l^) this is a 
substantial error. 

We proceed then to find a corrected version of Horton's 
law of stream numbers. Returning to equation (|ll]), we 
see that it is only valid in the limit Q — > oo. Defining 
n'(uj, Q) as the actual number of streams of order to within 
a basin of order f2, we have 

n' (w, ft)=2n'(uj + l,Q.) + ^2T v n'(u + v,Q,). (20) 

i/=i 
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This equation may be exactly solved. Considering the 
above expression for n'(u>, fi) and the corresponding one 
for n'(uj + 1, f2) we can reduce this to a simple difference 
equation, 

n(ui,n) = {2+Re+T 1 )n , (w+l,Q)-2R e n'(^+2,fl) (21) 

which has solutions of the form /i fc . Applying the con- 
straints that n'(Cl, Q) = 1 and n'(ft - 1, 0) = T x + 2, we 
obtain 



n'( W ^)= C ( M+ ) f2 - + (l- C )(/i_) S2 - 

where 

2fi ± = (2 + R e + Ti) ± [(2 + Rt + TO 2 - 8R e ] 



1/2 



and 



c = R n (R n - R e )/(R 2 n -2Re). 



(22) 
(23) 
(24) 



Note that R n — fi + and we will use the notation R n 
in place of This observation regarding Horton's law 
of stream numbers was first made by Tokunaga and 
later by Smart ]60| . In particular, Tokunaga noted that 
this would explain the deviation of Horton's law for the 
highest orders of a basin, a strong motivation for his 
work. 

We can now define an effective Horton ratio, R n '(oj, f2) 
as follows: 

Rn'(u},n) = n'(u-l,Si)/ri(u,n) 

= R n (l + 0(R;jR.J n ^A (25) 



The typical values of Horton's ratios in (19) give i?* = 1. 
In this case, R n {uj, f2) converges rapidly to R n with an 
error of around one per cent for u> = — 3. 



B. Hack's law 



As discussed in section [VC, an estimate of total 



drainage area of a basin is given by the total length of 
all streams within the basin. Summing over all stream 
orders and using the numbers n'(u>,Q,) given by equa- 
tions (0^) and (E3f) we have that 



* \S1 



ci(Rn) n + c 2 (Ri) u - c 3 (R* n ) 



(26) 



where c\ = c/(R n — Re), c 3 = (1 — c)/(Re — i?*) and 
C2 = C3 — ci with c being given in equation ( p4| ) . Slightly 
more complicated is the estimate of a(u>, f2), the drainage 
area of a basin of order u> within a basin of order f2: 

UJ' = 1 

= l/n'(uj,n) [ Cl (R n f(l - (Re/RnT) 
+ c 3 (R e r(R* n f-"(l - (K/RtT)] -(27) 



Now, for 1 <IC u> <C ^ (typically, 3 < ui < ft ~ 2 is 
sufficient), this expression is well approximated as 



a(w,fi) ~ (RnY 



(28) 



since R n > Re > R„- 

Thus, we have also shown here that R a = R n . While it 
is true that we would have obtained the same with a naive 
use of Horton's laws, we have both made the derivation 
thorough and established the correction terms found in 
equation (p7|) . This will be investigated further in the 
next section. 

Finally, usingthis result and the estimate oc (Ri)^ 
from equation (jq), it follows that 

L OC {ReY = (iky-In*/!**, „ (^lnJlf/lnft, (2g) 



which is precisely Hack's law. Comparing equations (|2£ 
and (||) , Hack's exponent is found in terms of the Horton 
ratios R n and Re as 



h = 



In Re 
In Rn 



(30) 



There is one minor caveat to the derivation in ( |29| ) and, 
for that matter, to most other derivations in this paper. 
Equation (^9|) only holds for the characteristic areas and 
lengths a u and l u . Since these quantities grow exponen- 
tially with u>, the derivation gives evenly spaced points 
on a log-log plot lying on a straight line. Clearly, this 
would indicate that the actual relationship is continuous 
and linear on a log-log plot. Indeed, there is no obvi- 
ous reason that a network would prefer certain lengths 
and areas. The averaging of stream lengths and areas 
brought about by the imposition of stream ordering nec- 
essarily removes all information contained in higher order 
statistics. Motivated by this observation, generalizations 
of the laws of Tokunaga, Horton and Hack to laws of 
distributions rather than averages is in progress |58[. 



VII. THERE ARE ONLY TWO HORTON 
RATIOS 

In deriving Hack's law in the previous section we ob- 
tained from equation ( pq) that R a = R n . This redun- 
dancy in Horton's laws is implicit in, amongst others, 
the works of Horton Q and Hack j5| but has never been 
stated outright. As noted previously, Peckham also ob- 
tains a similar result for a topological quantity, the num- 
ber of source streams in a basin, that is used as an es- 
timate of area. Thus, we see that for a landscape with 
uniform drainage density, Horton's laws are fully spec- 
ified by only two parameters R n and Re- This further 
supports our claim that Tokunaga's law and Horton's 
laws are equivalent since we have shown that there is an 
invertible transformation between (Ti,Rt), the param- 
eters of Tokunaga's law, and (R n ,Re) (equations JT^), 
( |l5j ) and (17)). In this section, we present data from real 
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Order u> for stream number data 
6 4 2 

O stream number 
□ stream area (km 2 " 1 




2 4 6 8 

Order to for stream area data 

FIG. 4: Average area and stream number as functions of 
stream order for Kentucky River, Kentucky (data taken from 
Peckham |^|). The stream number data is reversed for sim- 
pler comparison with the area data. The Horton ratios are 
estimated to be R n » R a — 4.65 ± 0.05. 



networks that support the finding R n = R a . We also 
address reported cases that do not conform to this re- 
sult and consider a possible explanation in light of the 
correction terms established in equation (26). 

Excellent agreement for the result R n — R a in real 
networks is to be found in the data of Peckham 1381 . The 



data is taken from an analysis of digital elevation models 
(DEM's) for the Kentucky River, Kentucky and the Pow- 
der River, Wyoming. Figure ^ shows average area and 
stream number plotted as a function of order for the Ken- 
tucky River while Figure ^] shows the same for the Powder 
river. Note that stream number has been plotted against 
decreasing stream order to make the comparison clear. 
The exponents R a and R n are indistinguishable in both 
cases. For the Kentucky river, R n « R a = 4.65 ± 0.05 
and for the Powder river, R n « R a = 4.55 ± 0.05. Also 
of note here is that the same equality is well satisfied by 
Scheidegger's model where numerical simulations yield 
values of R a = 5.20 ± 0.05 and R n = 5.20 ± 0.05. 

Note the slight deviation from a linear form for stream 
numbers for large lu in both cases. This upwards con- 
cavity is as predicted by the modified version of Horton's 
law of stream numbers for a single basin, equation (p2|). 

At the other extreme, the fit for both stream areas 
and stream numbers extends to lo = 1. While this may 
seem remarkable, it is conceivable that at the resolution 
of the DEM's used, some orders of smaller streams may 
have been removed by coarse-graining. Thus, lu = 1 may 
actually be, for example, a third order stream. Note that 
such a translation in the value of u> does not affect the 
determination of the ratios as it merely results in the 
change of an unimportant multiplicative constant. If lu t 



Order u for stream number data 

6 4 2 

O stream number 
□ stream area (km 2 " 1 




2 4 6 

Order lu for stream area data 



FIG. 5: Average area and stream number as functions of 
stream order for Powder River, Wyoming (data taken from 
Peckham H). Here the ratios are R n « R a = 4.55 ± 0.05. 



is the true order and lu = 
integer, then, for example, 



where m is some 



const x (R n ) u 



(31) 



This is only a rough argument as coarse-graining does 
not necessarily remove all streams of low orders. 

At odds with the result that R n = R a are past mea- 
surements that uniformly find R a > R n at a number of 
sites. For example, Rosso et al. in ^] examine eight river 
networks and find R a to be on average 40 % greater than 
R n . Clearly, this may be solely due to one or more of 
the our assumptions not being satisfied. The most likely 
would be that drainage density is not uniform. However, 
the limited size of the data sets points to a stronger pos- 
sibility which we now discuss. 

In the case of |5^] , the networks considered are all third 
or fourth order basins with one exception of a fifth order 
basin. As shown by equation (p6|), if Horton's laws of 
stream number and length are exactly followed for all or- 
ders, Horton's law of area is not obeyed for lower orders. 
Moreover, the former are most likely asymptotic relations 
themselves. It is thus unsatisfactory to make estimates of 
Horton's ratios from only three or four data points taken 
from the lowest order basins. Note that the Kentucky 
and Powder rivers are both eighth order networks and 
thus provide a sufficient range of data. 

We consider more precisely how the corrections to the 
scaling of area given in equation (^?j) would affect the 
measurement of the Horton ratios. Figure |^ shows an 
example of how stream number, length and area might 
vary with lu. It is assumed, for the sake of argument, that 
stream number and length scale exactly as per Horton's 
laws and that area behaves as in equation (E71) , satisfying 
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Order u> for stream number data 

10 8 6 4 2 



° stream number 
A stream area (km 2 ) 
n stream length (km) 




2 4 6 8 10 

Order u) for stream area and length data 

FIG. 6: An explanation for the empirical finding that R n < 
R a . Fitting a line to the stream area for only low uj would 
result in an overestimate of its asymptotic slope. For stream 
number, its slope would be underestimated. 



u) range 


1,2,3 1,2,3,4 1, 


2,3,4, 


5 4,5,6,7,8 


Rn 


2.92 3.21 


3.41 


3.99 


Ra 


5.29 4.90 


4.67 


4.00 



TABLE III: Values of Horton ratios obtained from least 
squares estimates of slopes for data represented in Figure ^| 
The range indicates the data points used in the estimate of the 
slopes. The ratios obtained from the low order data demon- 
strate substantial error whereas those obtained from the mid- 
dle data essentially give the true values of R n = R a = 4. 



Horton's law of area only for higher values of w. The plot 
is made for the example values R„ — 4 and Ri — 2. The 
prefactors are chosen arbitrarily so the ordinate is of no 
real significance. 

A measurement of R a from a few data points in the low 
lo range will overestimate its asymptotic value as will a 
similar measurement of R n underestimate its true value. 
Estimates of R n and R a from a simple least squares fit 



for various ranges of data are provided in Tabic II] 



Thus, the validity of the methods and results from past 
work are cast in some doubt. A reexamination of data 
which has yielded R a 3> R n appears warranted with an 
added focus on drainage density. Moreover, it is clear 
that networks of a much higher order must be studied to 
produce any reasonable results. 



VIII. FRACTAL DIMENSIONS OF 
NETWORKS: A REVISION 

A number of papers and works over the past decade 
have analyzed the relationships that exist between Hor- 
ton's laws and two fractal dimensions used to describe 
river networks §, §, f|, ||, ||, f| . These are D, the 
dimension which describes the scaling of the total mass 
of a network, and d, the fractal dimension of individual 
streams that comprises one of our assumptions. In this 
section, we briefly review these results and point out sev- 
eral inconsistencies. We then provide a revision that fits 
within the context of our assumptions. 

Our starting point is the work of La Barbera and 
Rosso j?]] which was improved by Tarboton et al. to 
give H 



D 



^Rn 

In Rt 



(32) 



We find this relation to be correct but that the assump- 
tions and derivations involved need to be redressed. To 
see this, note that equation (|32| ) was shown to follow from 
two observations. The first was the estimation of iV(^i), 
the number of boxes of size l\ x i\ required to cover the 
network M: 



In R n /\nR. e 



(33) 



where l\ is the mean length of first order stream seg- 
ments. Note that Horton's laws were directly used in 
this derivation rather than the correctly modified law of 
stream numbers for single basins (equation (p2|)). Nev- 
ertheless, the results are the same asymptotically. The 
next was the inclusion of our second assumption, that 
single channels are self-affine [||. Thus, it was claimed, 
l\ ~ S~ d where S is now the length of the measuring 
stick. Substitution of this into equation (p3|) gave 



N(S) ~ s- dlnRn / URe , 



(34) 



yielding the stated expression for D, equation (|32|). 

However, there is one major assumption in this work 
that needs to be more carefully examined. The network is 
assumed to be of infinite order, i.e., one can keep finding 
smaller and smaller streams. As we have stated, there is a 
finite limit to the extension of any real network. The pos- 
sible practical effects of this are pictorially represented 
in Figure ^. Consider that the network in question is 
of actual order f2. Then there are three possible scal- 
ing regimes. Firstly, for a ruler of length S ^> £i, only 
the network structure may be detected, given that indi- 
vidual streams are almost one-dimensional. Here, the 
scaling exponent will be \nR n /\nRe. Next, as S de- 
creases, the fractal structure of individual streams may 
come into play and the exponent would approach that 
of equation (|34|). Depending on the given network, this 
middle section may not even be present or, if so, perhaps 
only as a small deviation as depicted. Finally, the contri- 
bution due to the overall network structure must vanish 
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al. fB9 found that 
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FIG. 7: A schematic representing the problems associated 
with measuring the fractal dimension of a single river net- 
work. Here, the box counting method is assumed and 8, which 
has the units of length, is the side length of the N(5) boxes 
needed to cover the network. For box sizes much greater than 
i\ x li, only the network structure is detected while for box 
sizes smaller than £\x£\, the measurement picks out the frac- 
tal dimension of individual streams. Some deviation towards 
the scaling suggested by equation ( fil] ) may occur between 
these two limits. 



by the time S falls below l\ . From this point on, the mea- 
surement can only detect the fractal nature of individual 
streams and so the exponent must fall back to d. 

We therefore must rework this derivation of equa- 
tion As suggested in the definition of d in sec- 
tion [TVBl it is more reasonable to treat networks as grow- 
ing fractals. Indeed, since there is a finite limit to the 
extent of channelization of a landscape, there is a lower 
cut-off length scale beyond which most network quanti- 
ties have no meaning. The only reasonable way to exam- 
ine scaling behavior is to consider how these quantities 
change with increasing basin size. This in turn can only 
be done by comparing different basins of increasing order 
as opposed to examining one particular basin alone. 

With this in mind, the claim that equation ( |32] ) is the 
correct scaling can be argued as follows. Within some 
basin of order fi, take a sub-basin of order u>. Consider 
N(u), the number of boxes of side length l\ required 
to cover the sub-network. This is essentially given by 
the total length of all the streams in the network. This 
is given by the approximation of equation ( |2§| ) and so 
we have that N(u>) cx (i?„) w . Using the fact that £ UJ = 
{Ri^-Hx we then have that N(u) oc (4,/£i) lnfl " /lni? '<. 
The difference here is that l\ is fixed and pertains to the 
actual first order streams of the network. By assumption, 
we have that £ u oc L d and thus 



N(L) cx L dlaRn ^ Ul 



(35) 



which gives the same value for D as equation (p2[). 

There are two other relations involving fractal dimen- 
sions that also need to be reexamined. Firstly Rosso et 



J lnR a 



(36) 



Combining equations ( p2\ ) and (36), they then obtained 

. hi Rn 



D = 2- 



lnR a 



(37) 



However, equation (|36|) and hence equation fl37|) are both 
incorrect. 

There is a simple explanation for this discrepancy. In 
deriving equation ( |36| ) , Rosso et al. make the assumption 
that h = d/2, a hypothesis first suggested by Mandel- 
brot PJ. In arriving at the relation h = d/2, Mandelbrot 
states in jl| that "(basin area) 1 / 2 should be proportional 
to (distance from source to mouth as the crow flies)." 
In other words, a oc L 1 / 2 . However, as noted in equa- 
tion (^) , observations of real networks show that a oc L D 
where D < 2 Furthermore, on examining the re- 
sult h = In Ri I In R n with the expression for D in equa- 
tion (|33) we see that 



(38) 



which suggests that this hypothesis is valid only when 
D = 2. Consider also the test case of the Scheidegger 
model where h = 2/3, D — 3/2 and d = 1 (see Tabic ||). 
Using these values, we see that equation (B81) is exactly 
satisfied while the relation h — d/2 gives h = 1/2 ^ 2/3. 

Now, if h = d/D is used in place of h = d/2 in de- 
riving equation (|36| ) then equation (|32| ) is recovered. It 
also follows that equation ([57]) simplifies to the statement 
R a = R n , further demonstrating the consistency of our 
derivations. Thus, the two equations ( ^6| ) and ( p7| ) be- 
come redundant and the only connection between Hor- 
ton's ratios and network dimensions is given by equa- 
tion (||). 

An important point is that D < 2 does not imply 
that drainage basins are not space filling. This exponent 
shows how basin area changes when comparing different 
basins with different values of L, i.e., a oc L D . Any given 
single basin has of course a fractal dimension of 2. The 
equating of the way basin sizes change with the actual di- 
mension of any one particular basin is a confusion evident 
in the literature (see, for example, Q). Incorporating the 
effects of measuring basin area with boxes of side length 
5 in the relation a oc L D would lead to the form 



~2tD 



cll(5) oc 8 L 



(39) 



where the subscript L has been used to emphasize that 
different values of L correspond to different basins. Thus, 
for any given basin (i.e., for fixed L), the area scales with 
5 while for a fixed S, areas of different basins scale as per 
equation (||). 

It should also be emphasized that the relationship 
found here between Hack's exponent and the fractal di- 
mensions d and D is one that is explicitly derived from 
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the assumptions made. The observation that basin ar- 
eas scale non-trivially with L follows from these starting 
points and thus there is no need to assume it here. 



IX. OTHER SCALING LAWS 

We now address three remaining sets of scaling laws. 
These are probability distributions for areas and stream 
lengths, scaling of basin shape and Langbein's law. 

As introduced in equation (|J), probability distribu- 
tions for a and I are observed to be power law with ex- 
ponents r and 7 |^3| . Both of these laws have previousl 
been derived from Horton's laws. De Vries et al 
found a relationship between r, R n and Rg but did not 
include d in their calculations while Tarboton et al. || 
obtained a result for 7 that did incorporate d. 

Again, both of these derivations use Horton's laws di- 
rectly rather than the modified version of equation (22). 
Asymptotically, the same results are obtained from both 
approaches, 



Hence, 



lni^ 
In R„ 



and 



7 



In Rn 
In R f 



(40) 



Using the form of the Hack exponent found in equa- 
tion (pq ) and equation j32|), further connections between 
these exponents are found: 



t = 2 — h and 7 = — . 

h 



(41) 



One important outcome concerns the fact that only one 
of the exponents of the triplet (h,T, 7) is independent. 
Previously, for the particular case of directed networks, 
this has been shown by Meakin et al. |l4|] and further 
developed by Colaiori et al. [Q. Directed networks are 
those networks in which all flow has a non-zero positive 
component in a given direction. In a different setting, 
Cieplak et al. also arrive at this same conclusion for what 
they deem to be the separate cases of self-similar and self- 
affine networks although their assumptions arc that d < 1 
and D < 2 are mutually exclusive contrary to empirical 
observations J65|. In the case of non-directed networks, 
Maritan et al. have found one scaling relation for these 
three exponents, 7 = l+(r— l)/h and, therefore, that two 
of these three exponents are independent. They further 
noted that r = 2 — h is an "intriguing result" suggested 
by real data ||. In the present context, we have obtained 
this reduction of description in a very general way with, 
in particular, no assumption regarding the directedness 
of the networks. 

The scaling of basin shapes has been addressed already 
but it remains to show how it simply follows from our as- 
sumptions and how the relevant exponents are related. 
It is enough to show that this scaling follows from Hack's 
law. Now, the area of a basin is related to the longi- 
tudinal length L and the width L± by a — L±L, while 
the main stream length scales by assumption like I ~ L d . 



L d ~ (L ± L) h 



(42) 



where the fact that h = d/D has been used. Comparing 
this to equation (^) we obtain the scaling relation 



H = D-1. 



(43) 



The last set of exponents we discuss are those relating 
to Langbein's law ||. Langbein found that A, the sum 
of the distances (along streams) from stream junctions to 
the outlet of a basin, scales with the area of the basin. 
Recently, Maritan ct al. || introduced the quantity A, 
which is an average of Langbein's A except now the sum 
is taken over all points of the network. Citing the case 
of self-organized critical networks, they made the claim 
that 



AocL ¥ 



(44) 



Further, they assumed that ip — d although it was noted 
that there is no clear reason why this may be so since 
there are evident differences in definition (A involves dis- 
tances downstream while d involves distances upstream). 
We find this scaling relation to hold in the present frame- 
work. We further consider the two related quantities A 
and A, respectively the sum over all points and the av- 
erage over all junctions of distances along streams to the 
basin outlet. 

The calculations are straightforward and follow the 
manner of previous sections. We first calculate A(w, f2), 
the typical distance to the outlet from a stream of order 
w in an order £1 basin. Langbein's A, for example, is then 
obtained as Y^=i n(uj,Q)X(uj,Q). We find the same scal- 
ing behavior regardless of whether sums are taken over 
all points or all junctions. Specifically we find 

A- A~ a 1+lniWlnR ™ and A ~ A ~ L d (45) 

yielding the scaling relations 

/3 = P = 1 + ln^/ln-R„ and (p = ip = d. (46) 

Note that the second pair of scaling relations admit other 
methods of measuring d. The large amount of averaging 
inherent in the definition of the quantity A would suggest 
that it is a more robust method for measuring d than one 
based on measurements of the sole main stream of the 
basins. 

Maritan et al. || provide a list of real world measure- 
ments for various exponents upon which several com- 
ments should be made. Of particular note is the rela- 
tionship between r = 2 — h. This is well met by the 
cited values 1.41 < r < 1.45 and 0.57 < h < 0.60. Also 
reasonable is the estimate of h given by d/D (D = <f> in 
their notation) which is 0.58 < h < 0.65. 

The values of 7 and ip, however, do not work quite 
so well. The latter does not match d within error 
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bars, although they are close in absolute value with 
ip = 1.05 ± 0.01 and d = 1.10 ± 0.01. The length distri- 
bution exponent 7 may be found via 3 separate routes: 
7 = 1/h = D/d = 1/(2 — t). The second and third equal- 
ities have been noted to be well satisfied and so any one 
of the 3 estimates of 7 may be used. Take, for exam- 
ple, the range 0.58 < h < 0.59, which falls within that 
given by /i = 2 — r, h = d/D and the range given for h 
itself. This points to the possibility that the measured 
range 1.8 < 7 < 1.9 is too high, since using 7 = 1/h 
yields 7 = 1.74 ± .02. Also of note is that Maritan et 
al.'s own scaling relation y = l + (r — l)/h would suggest 
7 = 1.74 ±.05. 

Better genera l ag reement with the scaling relations is 
to be found in [p0[ in which Rigon et al. detail specific 
values of h, T and 7 for some thirteen river networks. 
Here, the relations r = 2 — h and j — 1/h are both well 
satisfied. Comparisons for this set of data show that, on 
average and given the cited values of h, both r and 7 are 
overestimated by only 2 per cent. 

X. CONCLUDING REMARKS 

We have demonstrated that the various laws, expo- 
nents and parameters found in the description of river 
networks follow from a few simple assumptions. Further, 
all quantities are expressible in terms of two fundamen- 
tal numbers. These are a ratio of logarithms of Horton's 
ratios, lni?„/lni?^, and the fractal dimension of individ- 
ual streams, d. There are only two independent param- 
eters in network scaling laws. These Horton ratios were 
shown to be equivalent to Tokunaga's law in informa- 
tional content with the attendant assumption of uniform 
drainage density. Further support for this observation is 
that both the Horton and Tokunaga descriptions depend 
on two parameters each and an invertiblc transformation 
between them exists (see equations ( |l~5| ) and (|l7j)). 
A summary of the connections found between the various 
exponents is presented in Table 

It should be emphasized that the importance of laws 
like that of Tokunaga and Horton in the description of 
networks is that they provide explicit structural informa- 
tion. Other measurements such as the power law proba- 
bility distributions for length and area provide little in- 
formation about how a network fits together. Indeed, 
information is lost in the derivations as the Horton ra- 
tios cannot be recovered from knowledge of lni?„/lni?£ 
and d only. 

The basic assumptions of this work need to be criti- 
cally examined. Determining how often they hold and 
why they hold will follow through to a greater under- 
standing of all river network laws. One vital part of any 
river network theory that is lacking here is the inclusion 



of the effects of relief, the third dimension. Another is 
the dynamics of network growth: why do mature river 
networks exhibit a self-similarity that gives rise to these 
scaling laws with these particular values of exponents? 
law: parameter in terms of R n , Re and d: 

T v — Ti(Rt) v 1 T\ = R n — Ri — 2 + IRej R n 
Rt = Re 

I ~ L d 
n w+ i/n u = Rn — 
L+i/L = Re 

lu+l/lu — Ri — 







Ra : 


— Rn 


I ~ 


a h 


h = 




a ~ 


L D 


D = 


= din Rn/ In R e 


L± <- 


•> L H 


H = 


-- dlnRn/lnR e - 1 


P(a). 


- a~ T 


r = 


2 - XnRn/lnRn 


P(l). 


- r 1 


7 


In Rn/ In R e 


A - 




'■ 


1 + ln^/ln Rn 


A ~ 




ip = 


d 


A- 




P = 


1 + ln R e /\nRn 


A ~ 




(p = 


d 



TABLE IV: Summary of scaling laws and the scaling relations 
found between the various exponents. Compare with table M. 



Also, extensive studies of variations in drainage density 
are required. The assumption of its uniformity plays a 
critical role in the derivations and needs to be reexam- 
ined. Lastly, in those cases where these assumptions are 
valid, the scaling relations gathered here provide a pow- 
erful method of cross-checking measurements. 

Finally, we note that work of a similar nature has re- 
cently been applied to biological networks The as- 
sumption analogous to network self-similarity used in the 
biological setting is considerably weaker as it requires 
only that the network is a hierarchy. A principle of mini- 
mal work is then claimed to constrain this hierarchy to be 
self-similar. It is conceivable that a similar approach may 
be found in river networks. However, a generalization of 
the concept of a hierarchy and perhaps stream ordering 
needs to be developed since a 'Tokunagic network' is not 
itself a simple hierarchy. 



Acknowledgements 

We are grateful to R. Pastor-Satorras, J. Pelletier, G. 
West, J. Weitz and K. Whipple for useful discussions. 
The work was supported in part by NSF grant EAR- 
9706220. 



[1] B. B. Mandelbrot, The Fractal Geometry of Nature (Free- [2] R. E. Horton, Bull. Geol. Soc. Am 56(3), 275 (1945). 
man, San Francisco, 1983). 



14 



[3] W. B. Langbein, U.S. Geol. Surv. Water-Supply Pap. W [35 
0968-C, 125 (1947). 

A. N. Strahler, Bull. Geol. Soc. Am 63, 1117 (1952). [36 
J. T. Hack, U.S. Geol. Surv. Prof. Pap. 294-B, 45 (1957). 
D. G. Tarboton, R. L. Bras, and I. Rodrfguez-Iturbe, [37 
Water Resour. Res. 24(8), 1317 (1988). 

P. La Barbera and R. Rosso, Water Resour. Res. 25(4), [38 
735 (1989). [39 

D. G. Tarboton, R. L. Bras, and I. Rodrfguez-Iturbe, 
Water Resour. Res. 26(9), 2243 (1990). [40 
A. Maritan, A. Rinaldo, R. Rigon, A. Giacometti, and [41 
I. Rodrfguez-Iturbe, Phys. Rev. E 53(2), 1510 (1996). [42 
L. B. Leopold and W. B. Langbein, U.S. Geol. Surv. Prof. [43 
Pap. 500-A, 1 (1962). 

R. L. Shreve, J. Geol. 74, 17 (1966). [44 
A. D. Howard, Geogr. Anal. 3, 29 (1971). 

C. P. Stark, Nature 352, 405 (1991). [45 
P. Meakin, J. Feder, and T. J0ssang, Physica A 176, 409 

(1991) . [46 
G. Willgoose, R. L. Bras, and I. Rodrfguez-Iturbe, Water 
Resour. Res. 27(7), 1685 (1991). [47 
G. Willgoose, R. L. Bras, and I. Rodrfguez-Iturbe, Earth 
Surf. Proc. Landforms 16, 237 (1991). [48 
G. Willgoose, R. L. Bras, and I. Rodrfguez-Iturbe, Water [49 
Resour. Res. 27(7), 1671 (1991). [50 
S. Kramer and M. Marder, Phys. Rev. Lett. 68(2), 205 

(1992) . [51 
R. L. Leheny and S. R. Nagel, Phys. Rev. Lett. 71(9), 
1470 (1993). [52 
T. Sun, P. Meakin, and T. J0ssang, Phys. Rev. E 49(6), [53 
4865 (1994). 

T. Sun, P. Meakin, and T. J0ssang, Water Resour. Res. [54 
30(9), 2599 (1994). [55 

E. Somfai and L. M. Sander, Phys. Rev. E 56(1), R5 
(1997). [56 
I. Rodrfguez-Iturbe and A. Rinaldo, Fractal River 
Basins: Chance and Self-Organization (Cambridge Uni- [57 
versity Press, Great Britain, 1997). [58 
P. Bak, How Nature Works: the Science of Self- Organized 
Criticality ( Springer- Verlag, New York, 1996). 
A. E. Scheidegger, International Association of Scientific [59 
Hydrology Bulletin 12(1), 15 (1967). 

A. E. Scheidegger, Theoretical Geomorphology (Springer- [60 
Verlag, New York, 1991), third ed. [61 
A. N. Strahler, EOS Trans. AGU 38(6), 913 (1957). [62 
M. A. Melton, J. Geol. 67, 345 (1959). 

D. M. Gray, J. Geophys. Res. 66(4), 1215 (1961). [63 
R. Rigon, I. Rodrfguez-Iturbe, A. Maritan, A. Gia- 
cometti, D. G. Tarboton, and A. Rinaldo, Water Resour. 
Res. 32(11), 3367 (1996). [64 
J. E. Mueller, Geological Society of America Bulletin 83, 
3471 (1972). [65 
M. P. Mosley and R. S. Parker, Geological Society of 
America Bulletin 84, 3123 (1973). 

J. E. Mueller, Geological Society of America Bulletin 84, [66 
3127 (1973). 

P. La Barbera and R. Rosso, Water Resour. Res. 26(9), 
2245 (1990). 



E. Tokunaga, Geophys. Bull. Hokkaido Univ. 15, 1 
(1966). 

E. Tokunaga, Geogr. Rep., Tokyo Metrop. Univ. 13, 1 
(1978). 

E. Tokunaga, Trans. Jpn. Geomorphol. Union 5(2), 71 
(1984). 

S. D. Peckham, Water Resour. Res. 31(4), 1023 (1995). 
W. I. Newman, D. L. Turcotte, and A. M. Gabrielov, 
Fractals 5(4), 603 (1997). 

A. E. Scheidegger, Water Resour. Res. 4(5), 1015 (1968). 

J. W. Kirchner, Geology 21, 591 (1993). 

S. A. Schumm, Bull. Geol. Soc. Am 67, 597 (1956). 

H. Takayasu, I. Nishikawa, and H. Tasaki, Phys. Rev. A 

37(8), 3110 (1988). 

M. Takayasu and H. Takayasu, Phys. Rev. A 39(8), 4345 
(1989). 

H. Takayasu, Physcial Review Letters 63(23), 2563 
(1989). 

H. Takayasu, Fractals in the Physical Sciences (Manch- 
ester University Press, Manchester, 1990). 
H. Takayasu, M. Takayasu, A. Provata, and G. Huber, J. 
Stat. Phys. 65(3/4), 725 (1991). 

G. Huber, Physica A 170, 463 (1991). 

A. D. Abrahams, Water Resour. Res. 20(2), 161 (1984). 
D. R. Montgomery and W. E. Dietrich, Science 255, 826 
(1992). 

D. G. Tarboton, R. L. Bras, and I. Rodrfguez-Iturbe, 
Water Resour. Res. 25(9), 2037 (1989). 
R. L. Shreve, J. Geol. 75, 178 (1967). 
P. Haggett and R. J. Chorley, Network Analysis in Ge- 
ography (Edward Arnold, London, 1969). 
V. Gardiner, Geogr. Ann. 55A(3-4), 147 (1973). 
M. E. Morisawa, Geological Society of America Bulletin 
73, 1025 (1962). 

H. de Vries, T. Becker, and B. Eckhardt, Water Resour. 
Res. 30(12), 3541 (1994). 

W. S. Clock, The Geogr. Rev. 21, 475 (1931). 

P. S. Dodds, R. Pastor-Satorras, and D. H. Rothman, 

Fluctuation in River network scaling laws (1999), in 

preparation. 

R. Rosso, B. Bacchi, and P. La Barbera, Water Resour. 
Res. 27(3), 381 (1991). 

J. S. Smart, Water Resour. Res. 3(3), 773 (1967). 

J. Feder, Fractals (Plenum Press, New York, 1988). 

P. La Barbera and G. Roth, Hydrol. Processes 8, 125 

(1994). 

C. P. Stark, Stream networks on a Bethe lattice: Cayley 
trees, invasion percolation and branching ratios (1997), 
preprint. 

F. Colaiori, A. Flammini, A. Maritan, and J. R. Banavar, 
Phys. Rev. E 55(2), 1298 (1997). 

M. Cieplak, A. Giacometti, A. Maritan, A. Rinaldo, 

I. Rodrfguez-Iturbe, and J. R. Banavar, J. Stat. Phys. 
91(1/2), 1 (1998). 

G. B. West, J. H. Brown, and B. J. Enquist, Science 276, 
122 (1997). 



